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We present a new mechanism for air shower radio emission due to the sudden absorption of 
secondary particles when the shower front hits the ground. The electrons present in excess during 
the air shower development imply a net residual negative charge in the shower front. We show that 
for showers hitting the ground before the complete extinction of their electromagnetic component, 
the sudden vanishing of the net residual negative charge generates an electric field contribution 
in the kHz — MHz range. We characterize this radio contribution as a function of primary energy, 
arrival direction and antenna position, using the simulation code SELFAS2. We discuss the interest 
of this new predicted signal on detection and analysis of ultra-high energy cosmic-rays and we argue 
that the region in the shower of maximum emission of the electric field should not coincide with the 
region of maximum development. 
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I. INTRODUCTION 



through: 



The electromagnetic particles of air showers reaching 
the ground are detected by ground particles detectors 
such as Water Cerenkov Detector (used at the Pierre 
Auger Observatory for instance [T) or plastic scintil- 
lators (used in the Telescope Array experiment for in- 
stance p]). The measurement of their ground density 
allows an estimation of the core position and the pri- 
mary energy. These particles should be also detectable 
through the electric field emitted consecutively to their 
sudden absorption in the ground. We propose in this 
paper a complete characterization of this signal using 
simulations. We use the code SELFAS2 [3 , where each 
secondary electron and positron of the shower front is 
considered as a moving source. The total electric field 
emitted by the complete shower at any observation point 
x, is obtained after superposition of all individual contri- 
butions. Considering only the radiative part of the field 
(we neglect the first term of Eq. (17) in [3]), the total 
field detected by an observer located at x and at time t 
is given by: 
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where rj is the air refractive index, and B4 are the line 
of sight and the distance between the observer and the 
particle i, (Hi the velocity of this particle and its electric 
charge. The summation is done over the total number N t 
of particles that emitted an electric field detected by the 
observer at time t. All these quantities are evaluated at 
the retarded time t ret , related to the observer's time t 
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The Earth's magnetic field induces a systematic opposite 
drift of electrons and positrons during the shower devel- 
opment; this generates a residual current j± = qc(3± 
perpendicular to the shower axis. The variation of this 
current with the number of particles creates the main 
electric field contribution, in the MHz range. In addi- 
tion, the charge excess variation due to positrons annihi- 
lation and knock-on electrons during the air shower de- 
velopment (Askaryan effect in air [1]) gives a secondary 
contribution to the total electric field, also in the MHz 
range (see [3j for a complete description). Experimen- 
tal evidence for this specific contribution can be found 
in [51 [6] . Finally, the air refractive index alters the radio 
signal for observers located close to the shower axis (less 
than 200 m for a vertical shower, see [7])- The resulting 
electric field generated during the shower development in 
the atmosphere creates the principal pulse (PP), as al- 
ready observed since years by several experiments such as 
CODALEMA 019], LOPES QUE] or AERA [12 for in- 
stance. Since the last decade, the convergence of various 
experimental results and various theoretical approaches 
shows that the mechanism of the radio signal emitted by 
the shower during its development in the atmosphere is 
quite well understood. 

We suggest a new mechanism generating a clear radio 
signal due to the absorption of the shower front when 
reaching the ground. In this paper, we first describe the 
mechanisms at the source of this electric field. Then, we 
present the characteristics of this signal using simulated 
showers with SELFAS2, in both time and frequency do- 
mains. The dependence with the arrival direction, energy 
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of the primary cosmic ray and the relative position of the 
observer with respect to the shower core is discussed. Fi- 
nally, we discuss the importance of this new predicted 
signal in the field of ultra-high energy cosmic rays. 



II. SUDDEN DEATH RADIO EMISSION 

Depending on the air shower characteristics and on the 
altitude of the observation site, the number of secondary 
particles reaching the ground N(X gI d) can vary strongly. 
For a given ground altitude, A(A gr( j) mainly depends on: 

• the energy of the primary cosmic ray inducing the 
cascade in the atmosphere; 

• the air shower arrival direction (zenith angle); 

• the first interaction length X\ of the primary cos- 
mic ray; 

• the nature of the primary cosmic ray. 

Fig. [I] (top) shows the number of particles reaching the 
ground as a function of arrival direction, for showers in- 
duced by proton and iron nuclei of various energies. The 
corresponding longitudinal profiles are presented in Fig. [I] 
(bottom). The maximum number of particles iV max in 
the shower is obtained, by definition, at the depth of 
maximum development X max . We see that for a primary 
energy above 10 18 eV (1 EeV) and for arrival directions 
9 ^ 50°, iV(Xg r d) is not negligible with respect to iV max 
(A(A grd )/A max > 10%). 

A large majority of secondary electrons and positrons 
(at least 97%) for X ^ A max , have an energy smaller 
than 1 GeV (see [17]). Consequently, when the shower 
front hits the ground, the secondary electrons and 
positrons are stopped on a distance smaller than 1 m at 
1 GeV in a silicon medium |18j , similiar to the ground at 
the Auger site. The radiation of electrons and positrons 
under ground level is strongly attenuated by the small 
conductivity of this medium: cr^^nd = 8.10~ 3 S/m at 
80 MHz [T5], implying a complete attenuation of the ra- 
diation after a few centimeters. 

The very fast absorption of these particles — we call 
this phenomenon the sudden death (SD) of the shower — 
is at the origin of the emission of an electric field through 
various mechanisms: 

• the sudden disappearance of the net negative 
charge in the shower front generated by secondary 
electrons in excess, and the corresponding longitu- 
dinal current; 

• the sudden disappearance of the transverse current 
due to systematic opposite drift of electrons and 
positrons in the geomagnetic field; 
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FIG. 1: Top: number of particles reaching the ground JV(X gr< j) 
as a function of zenith angle, for showers initiated by proton and 
iron nuclei of various energies, computed for an altitude of 1400 m 
corresponding to the Auger site. The Greisen-Iljina-Linsley param- 
eterization 1 1 3 1 114| of the longitudinal profile is used to estimate 
7V(X grt j). The different lines (solid, dotted, dashed...) represent 
the number of ground particles, obtained for showers with a value 
of Xi set to the mean value < X\ (E, A) > according to the pri- 
mary energy and mass (using the QGSJET cross section extrap- 
olations 1 151 116| ). The fluctuation of Xi implies a fluctuation on 
7V(X g rd), represented by the variable color intensity. More intense 
means more probable. Bottom: corresponding longitudinal pro- 
files. The vertical gray lines represent the depth crossed by the 
shower until it reaches the ground, depending on the value of the 
zenith angle 8. 



• the asymetry in the ground distribution of electrons 
and positrons, induced by the charge separation. 

This electric field is emitted by a region close to the 
shower core. An observer will detect the signal at a 
time ~ d/c after the impact on the ground, where d is 
the distance to the shower core. The electric field emitted 
during the development of the shower in the atmosphere 
creates the PP as discussed in the introduction, which is 
detected before this sudden death pulse (SDP). The time 
interval between the PP and the SDP depends on the ge- 
ometry (the air shower arrival direction and the antenna 
position with respect to the ground shower core) but also 
on the instant of the maximum of emission during the air 
shower development. 
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In the next section, we quantify and study in details 
the SDP using the simulation. 



III. DESCRIPTION OF THE PULSE 
Example of vertical simulated showers 

We ran simulations of proton-induced vertical showers 
falling at the origin of the coordinate system for energies 
0.1, 0.3, 1, 3, 10, 30, 100 EeV with the code SELFAS2 for 
48 antennas located between 100 m and 800 m by steps of 
50 m at the geographic east, west and north of the shower 
core and one antenna located at the origin. The ground 
level is set at 1400 m and we consider the geomagnetic 
field measured at the Pierre Auger Observatory. The first 
interaction point (X\) at a given energy is extracted from 
the QGSJET data as in Fig.[l] and is taken as a function 
of energy as 52.8,50.8,48.6,46.8,44.9,43.2,41.5 g cm" 2 , 
respectively. The electric field received by each antenna 
as a function of time is calculated in the EW, north-south 
(NS) and vertical (V) polarizations. The observer's origin 
of time t = corresponds to the time when the shower 
axis hits the ground at the core position (0, 0, 1400). 

Figs. [2] and [3] presents the EW and V electric fields, 
respectively, as a function of time at 500 m and 600 m at 
the east of the shower core. 
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FIG. 3: Same as Fig.[2]for the V polarization. 



cal plane containing the shower axis and the observa- 
tion point, in agreement with an electric field aligned 
with n — j3 (see Eq. [I]). For a vertical shower, an ob- 
server measuring only the electric field in the horizontal 
plane will reconstruct a radial polarization pattern, in 
agreement with Eq. [I] which predicts a radial polariza- 
tion pattern in the horizontal plane when the shower is 
vertical ((3 ~ — z). We checked that the contribution to 
the SDP is dominated by the excess of negative charge 
in the shower. 



> 
E 



0.05 
-0.00 



| -0.05 
I -0.10- 

CD 

0.15 



axis distance 500 m 
axis distance 600 m 
3.00 EeV 
10.0 EeV 



500 1000 1500 2000 
time (ns) 



2500 3000 



FIG. 2: EW electric field vs time for two antennas located at 500 m 
and 600 m at the east of the shower core, for vertical showers at 3 
and 10 EeV. 



The PP due to the various mechanisms generating the 
electric field in the shower during the development in the 
atmosphere is clearly visible at the beginning of the trace, 
with a maximum amplitude occuring at a time close to 
100 ns in this example. The SDP appears at roughly 
1600 ns and 2000 ns for antennas located at 500 m and 
600 m of the shower core, respectively. Note that in this 
example, no SDP is visible in the NS polarization for 
these antennas located at the east of the shower axis. 
This means that the SD electric field lies in the verti- 



Time structure of the SDP 

From the simulation, the time interval between the 
time origin and the time tg^p when the SDP reaches its 
maximum value is simply given by St = d/c, in agreement 
with the hypothesis that the SDP is due to the disap- 
pearance of secondary particles when hitting the ground. 
This result is also verified for all antennas used in this 
simulation. 

We can infer the shape of the SDP. For that, we de- 
scribe the shower front as a macroscopic spatial charge 
density defined by Q(x, y, z, t rct ). This function describes 
the thickness and lateral spread of the shower front. Be- 
fore reaching the ground, the time variation of Q due to 
the evolution of the number of particles in the shower is 
small compared to the charge variation induced by the 
sudden absorption of the particles in the ground, which is 
a quasi instantaneous event. The time derivative of a de- 
creasing charge to zero explains why the SDP is monopo- 
lar, as seen in Figs [2] and [3j When the shower front hits 
the ground, the observer sees the shower front disappear- 
ing at the speed of light. The duration of the complete 
absorption in the ground depends on the inclination of 
the shower, on its curvature and its thickness. The thick- 
ness of the shower front (« 2 m) is negligible compared 
to the size of the lateral extension (Moliere radius ~ 100 
m) , so that we expect that the time structure of the SDP 
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is mainly due to the lateral spread of the shower front. 
The SDP time structure provides an image of the parti- 
cle ground footprint. For this reason, the shape of the 
pulse should be predicted quite easily by simple geomet- 
rical considerations. Two equivalent points of view can 
be adopted to model the shape of the SDP. 



Analytical description using the LDF, for a vertical shower 

Let's first consider the case of a vertical shower where 
the ground particles density (per unit area) is well de- 
scribed by a symetric lateral distribution function (LDF) 
C{p), where p is the distance to the shower axis. The 
shower core is located at the position (x c ,y c ) and the 
observer at the position (x a ,y a ). We call d c the dis- 
tance between the shower core and the observer and a 
the angle under which the observer in A sees the shower 
core (see Figj4]). Let M be a point on the ground lo- 
cated at a distance r from A. At M, the particle den- 
sity is given by the value of the LDF at this point, at a 
distance p(r, 4>) — y/r 2 + d 2 — 2rd c cos(a — 4>) from the 
shower core. The number of particles around A in the 
interval [r; r + dr] is given by: 



AN 



(r, d c ) = r 



£(p(r,$)d0 
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According to Eq. [T] the particles located in the interval 




FIG. 4: System of coordinate used to estimate the time structure 
of the SDP. A is the observer position and M a point on the ground 
where the particle density is given by the LDF C(p), p which is the 
distance between M and the ground shower core. 

[r;r + dr] create an electric field -Esdp^a) observed by 
A proportional to 1/r, detected at the time tA — r/c: 



SDP 



(t A = r/c) oc 
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r dr 
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Using the Greisen lateral distribution function C(p) given 
in 20J at the ground level, we show the result obtained 
in Fig. [5] in the case of a vertical shower. The absolute 
amplitude of the analytic SDP model is fitted to the sim- 
ulated SDP obtained at 200 m and the same factor is 
used for the antenna located at 400 m from the shower 



core. The analytical model reproduces very well the time 
structure and the relative amplitude of the SDP. 



Description using the ground particles distribution 

The analytical model can be extended to inclined show- 
ers, provided we know the actual LDF, including the 
asymetries due to the shower inclination: the LDF de- 
pends on p but also on the relative position around the 
core. A more direct method to compute the expected 
SDP for inclined showers, is to use directly the infor- 
mation provided by the simulation. In the simulation, 
we record the time and ground position (ti,Xi,yi, Zi = 
1400 m) of each simulated particle i hitting the ground. 
We assume that these particles have the same velocity 
parallel to the shower axis (3. According to Eq.[IJ a given 
ground particle i located at the point M of coordinates 
(t.j, Xi,yi, Zi — 1400 m) creates an electric field at the po- 
sition A proportional to (n — f3)/Ri where Ri = \ \MA\\ 
and n = MA/R L . This signal is detected by the observer 
at the time tf = U + Ri/c. For this observer, the total 
signal detected in the time interval [tf;tf + St[ is the 
summation of all contributions of the ground particles 
verifying tf < ti + Ri/c < tf + St. This simple mod- 
elization of the SDP can be compared to the result of the 
simulation. In Fig. [5] we show the simulated electric field 
for two antennas located at 200 m and 400 m at the east 
of the shower axis for a vertical shower of 1 EeV. In Fig|6j 
we show the same analysis for an inclined shower coming 
from the north-east (6 = 30°, <f> = 45°) in the NS polar- 
ization for two antennas at the north of the shower core, 
located at 168 m and 475 m, respectively (axis distances 
157 m and 445 m). We see that the agreement is very 
good. For a given simulated shower (cither vertical or in- 
clined as shown in Figs [5] and [6]) , the absolute amplitude 
of the SDP model is fitted to the simulated SDP and the 
same factor is used for all antennas corresponding to this 
shower. This means that the model reproduces correctly 
both the time structure and the relative amplitude of this 
signal. 

For a given antenna, Fig. [7] shows that the model pre- 
dicts correctly the SDP in the three polarizations. The 
data presented here correspond to the simulated electric 
field computed at 268 m at the east of the shower core, 
for an inclined shower 9 = 30° coming from the north- 
east, at 1 EeV. The SDP model is computed for the three 
polarizations and the same normalization factor is used, 
this confirms that the SDP pulse is well described by an 
electric field along n (3. 



SDP amplitude characterization 

In this section we characterize the SDP amplitude A 
as a function of energy and core distance. We estimate 
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FIG. 5: SDP model applied to a vertical shower at 1 EeV using 
two antennas in the EW polarization at 200 m and 400 m at the 
east of the shower core. The red dots are the predicted SDP using 
the model based on the ground particles distribution provided by 
the simulation. The blue squares correspond to the SDP predicted 
using the analytical Greisen formula for the LDF. In both cases, 
the SDP is very well reproduced. 
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FIG. 6: SDP model applied to an inclined shower at 1 EeV coming 
from the north-east (9 = 30°, cj> = 45°) using two antennas in the 
NS polarization at 157 m and 445 m from the axis. The antennas 
are located at the east of the shower core, at 168 m and 475 m 
respectively. 



the amplitude by fitting the SDP with the function: 

U _ y-max I ■ 

r '-sdp I 



SDP(t) = .4cxp - 



(5) 



where tg^p is the time when the SDP reaches its maxi- 
mum value and r± is a time constant which can take dif- 
ferent values on the two intervals t ^ t'^p and t > t'^p. 



Vertical case 

Fig. [H] presents the dependence of A with the primary 
energy, for observers located at different core distances, 
at the east of the shower axis. The depence is almost 
linear. 



FIG. 7: SDP model applied to the three polarizations of a single 
antenna located at 268 m at the east of the shower core. The same 
normalization factor is used for the three polarizations. The shower 
is coming from 8 = 30°, <f> = 45°, the primary energy is 1 EeV. 
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FIG. 8: SDP amplitude in the EW polarization as a function of 
energy for different core distances. The amplitude varies almost 
linearly with the energy. 



Fig.[9]presents the dependence of A with core distance, 
at different energies. We see that A scales as the inverse 
of the core distance, proving the radiative nature of the 
SD emission as inferred in the previous section. We can 
finally propose a complete parameterization of the SDP 
amplitude A as a function of energy and core distance: 



100 m 



E 



1 EeV 



/iVm 



(6) 

the factor »4norm being the normalized amplitude in a 
given polarisation ("pol" being EW, NS or V) at 1 EeV 
and a reference core distance of 100 m. For instance, 
in the EW polarization for observers at the east of the 
shower axis, -A^orm = 15-6 /i(Vm -1 . We find the same or- 
der of amplitude for antennas at the west of the shower 
core in the EW polarization and at the north in the NS 
polarization. Note that the SDP amplitude is relatively 
higher in the V polarization than in the horizontal polar- 
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FIG. 9: SDP amplitude in the EW polarization as a function of 
core distance for different energies. The amplitude scales as 1/d. 



izations when compared to the PP amplitude. At fixed 
energy, the SDP amplitude decreases less slowly than 
the PP amplitude so that the ratio of the amplitudes 
SDP/PP can even be greater than 1 (but at energies 
above 20 EeV for horizontal polarizations). 



Inclined case 

At fixed energy, first interaction point and primary cos- 
mic ray, the number of particles reaching the ground is 
maximum when the shower is vertical (see Fig. [I]). There- 
fore, we expect a smaller amplitude for inclined showers. 
The amplitude scales in the same way as vertical showers 
(see Eq. [6|, in both energy and core distance; the refer- 
ence amplitude at 1 EeV for a core distance of 100 m is 
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of the order of 5 fiV/ui for a zenith angle of 30°. Fig 
presents the variation of the SDP amplitude A for an in- 
clined shower with = 30°, <f> = 45°, as a function of 
energy for different core distances. 



Frequency domain 



FIG. 10: SDP amplitude in the NS polarization as a function of 
energy for different core distances. The antennas are located at the 
north of the shower core. The SDP amplitude varies almost linearly 
with the energy. This case corresponds to an inclined shower of 
energy 1 EeV, with 9 = 30° and = 45° . 
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FIG. 11: PSD of the PP in the EW polarization for different 
core distances located at the east of the shower axis. This PSD 
corresponds to a vertical shower of energy 1 EeV. 



to a frequency of ~ 20 MHz whatever the core distance 



is, as can be seen in Fig. 12 



Contrarily to the time structure of the PP, that of the 
SDP does not depend on the primary energy nor on the 
distance to the shower core as can be seen in Figs. [2] 
and [3] Only the amplitude is affected by these parame- 
ters. Therefore, we expect the shape of the power spec- 
tral density (PSD) to be similar at all core distances and 
all energies as a function of frequency, the only difference 
being the normalization of the PSD. Fig. [IT] presents the 
PSD of the PP for a primary energy of 1 EeV and core 
distances varying from 100 m to 800 m. We observe, as 
already predicted by several simulation codes, that the 
PP PSD is coherent up to a certain frequency depending 
on the axis distance of the observer. This is not the case 
for the SDP since we observe that its PSD is coherent up 



DISCUSSION AND CONCLUSION 

We have characterized the electric field produced by 
the secondary e ± when they reach the ground. The am- 
plitude of this signal in the horizontal polarization, is 
of the order of 15 /LtV/m at 100 m of the shower core 
for a vertical shower, and a primary energy of 1 EeV. 
The amplitude scales quasi-linearly with the energy and 
decreases as 1/d, where d is the distance between the 
shower core and the observer. The SDP amplitude de- 
pends on the distance to the shower core, contrarily to 
the PP due to the secondary e during the shower de- 
velopment, which strongly depends on the distance to 
the shower axis. The polarization is oriented along the 
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FIG. 12: PSD of the SDP in the EW polarization for different 
core distances located at the east of the shower axis. This PSD 
corresponds to a vertical shower of energy 1 EeV. 



vector n — /3, where n is the normalized vector between 
the shower core and the observer position and /3 can be 
taken as the shower axis direction. The signal exists if 
the number of e in the shower at the ground level is 
high enough: for a given observation site, the condition 
to observe the SDP depends on the zenith angle and pri- 
mary energy of the shower. Fig. 13 presents the expected 



SDP amplitude at 100 m at the east of the shower core, 
in the EW polarization, as a function of zenith angle and 
energy. We assumed proton-initiated showers with the 
first interaction point X\ set to the average value ex- 
tracted from QGSJET data as a function of primary en- 
ergy. At the Pierre Auger Observatory site for instance, 
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FIG. 13: Expected amplitude of the SDP at 100 m at the east of 
the shower core, in the EW polarization, as a function of primary 
energy and zenith angle, computed for the site of the Pierre Auger 
Observatory. 



the ground altitude of 1400 m implies that a total of 
~ 10 9 e can reach the ground for a vertical shower of 



1 EeV. The SDP power spectral density is important at 
relatively low frequencies, below 20 MHz. If we are able 
to detect a signal at these frequencies, then it would be 
quite simple to trigger on actual cosmic rays because any 
trace having two pulses separated by some /j,s could be 
an excellent candidate. This would be a very specific sig- 
nature of cosmic rays and would help a lot in discarding 
background events. The SDP is expected to arrive at 
the observer location at a delayed time d/c with respect 
to the core time. The core is therefore located on a cir- 
cle centered on the observer position and of radius d/c. 
Using the information from several detectors, the shower 
core is at the intersection of the corresponding circles. 
The SDP amplitude is proportionnal to the total number 
of secondary e ± and also reflects their complete ground 
distribution (see Eq. [1]) , contrarily to particle detectors 
that sample at specific locations this ground distribution. 
Experimentally, one should be interested in recording the 
electric field over a duration greater than some /is to be 
able to observe the SDP (up to ~ 3 /j,s after the impact 
time at 1 km from the shower core). The antennas used 
should be also sensitive to horizontal directions. 

Because the SDP is generated by the end of the shower, 
its detection provides an absolute timing of the shower in 
the trace recorded by an antenna. For antennas not too 
close from the shower axis (because the effect of the air 
refractive index is important close to the shower axis), 
we expect a one-to-one correspondence [7] between the 
time in the antenna trace and the position of the source 
in the shower. This bijective relation can be established 
by simple geometrical considerations, if we assume that 
the shower front is point-like and moves at the speed 
of light c. Let tj be the time of the bin i in the time 
trace and ^sDp the time of the maximum amplitude of 
the SDP. The time f§DP corresponds to the end of the 
shower development, so that t t corresponds to the instant 
when the shower front was at a distance £(U) from the 
shower core, along the shower axis, given by: 



t(ti 



2 Atf - 2c Mi \\AC\ 



2(cAti 



\AC\\ + AC ■ u) 



(7) 



where Ati = igDp — U, AC is the vector from the antenna 
position A to the shower core C and u is the unit vector 
aligned with the shower axis (sin 8 cos <j>, sin sin cj>, cos ff) . 
If we consider the particular time bin corresponding to 
the maximum amplitude of the PP, U = tpp x , we can 
estimate the position — £(tpp x ) along the shower 

axis corresponding to the maximum production of elec- 
tric field, for this antenna. Wc finally compute the cor- 
responding altitude above sea level zgiax taking into ac- 
count the Earth's curvature, and we finally obtain the as- 
sociated depth X-j£°x, assuming the US standard model 
for the atmosphere [21] 



Fig. 14 shows the estimation 



of as a function of the distance to the shower 

core, for antennas located at the east, west and north 
of the core, for a vertical shower at 1 EeV. The error 



bars were computed by Monte-Carlo assuming an error 
on At = f§DP — ipp* varying linearly from 10 ns at 100 m 
of the shower core up to 20 ns at 1000 m of the shower 
core. In this figure, we see two important characteristics: 



1000 





800 - 


1 




g 


600 ' : 


&s 






400 ' : 




200 ' : 




o : 




200 400 600 800 1000 
core distance (m) 



^prod 
L max 

shower core, for antennas located at the east, west and north of 
the core. The red dashed line indicates the X max = 749 g cm -2 
for this vertical simulated shower at 1 EeV. The black dashed line, 
indicate the depth corresponding to the first inflexion point of the 
longitudinal profile, Xi n f = 525 g cm~ 2 . 

a convergence of the X^°£ estimation with increasing 
core distances and a systematic offset between the tree 
curves corresponding to antennas located at the north 
(blue), east (red) and west (green) of the shower core. 
The convergence can be explained by the fact that for an 
antenna close to the shower axis, the air refractive index 
generates a time dispersion of the signal which destroys 
the one-to-one correspondence between the time in the 
antenna trace and the corresponding source position in 
the shower (i rot is a non-monotonic function of t, see [7]). 
In that case, the relation given in Eq. [7J used to estimate 
IP™ d cannot be applied and the estimation of the 
is not valid. When the axis distance increases, the re- 
lation between t ret and t becomes monotonic and Eq. [7J 
can be used to estimate JP' a ° x d . The systematic offset 
between the three directions north, east and west is due 
the systematic difference in the position in the time trace 
of the PP. In this example of a vertical shower at 1 EeV 
with the geomagnetic field of the Pierre Auger Observa- 
tory site, the PP is systematically delayed for antennas 
located at the west of the shower axis with respect to an- 
tennas located at the north and east of the shower axis: 
^ppfwest > *pp X cast,„orth at a fixed core distance. The rea- 
son for that is the following: the two main mechanisms at 
the origin of the PP (respectively charge excess and trans- 
verse current variation during the shower development) 
produce two electric fields E n (aligned along n) and 

E e ± (aligned along (5\ ) that interfer constructively 
Pi 

or destructively following the antenna position with re- 



spect to the shower axis and to the geomagnetic field. In 
the configuration adopted in this paper, the interference 
is destructive at the west of the shower core, construc- 
tive at the east of the shower core, implying a systematic 
relative time shift between the time of maximum of the 
PP amplitudes between antennas located at the east and 
west of the shower core. This explains the systematic 
shift to higher values of JP r a ° x d for antennas at the west 

: clearly shows that the JP™ d is 
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of the shower core. Fig. 
not in agreement with the depth X max where the number 
of particles in the shower is maximum. The maximum 
of radio emission is close to the depth Aj n f at which the 
production rate of particles in the shower per g cm~ 2 is 
maximum (it corresponds to the first inflexion point of 
the longitudinal profile, [dN/dX] x ^ f = 0). 
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